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We construct a field theory to describe energy averaged quantum statistical properties of systems 
which are chaotic in their classical limit. An expression for the generating function of general 
statistical correlators is presented in the form of a functional supermatrix nonlinear cr-model where 
the effective action involves the evolution operator of the classical dynamics. Low-lying degrees 
of freedom of the field theory are shown to refiect the irreversible classical dynamics describing 
relaxation of phase space distributions. The validity of this approach is investigated over a wide 
range of energy scales. As well as recovering the universal long-time behavior characteristic of 
random matrix ensembles, this approach accounts correctly for the short-time limit yielding results 
which agree with the diagonal approximation of periodic orbit theory. 



I. INTRODUCTION 

The quantum description of systems which are chaotic in the classical limit is the subject of "quantum chaos" . 
A wide variety of physical systems fall into this category. Amongst those most commonly studied are the neutron 
resonances of atomic nuclei Q, Rydberg atoms in strong magnetic fields and electrons in weakly disordered 
metallic grains (quantum dots) The energy spectrum as a whole is specific for each individual chaotic system. 
However, in contrast to integrable systems, each eigenstate is characterized only by its energy rather than by a set 
of quantum numbers Q. The variables in the corresponding Schrodinger equation do not separate, and an analytical 
solution is prohibited. Therefore, a useful description of highly excited eigenstates of chaotic systems is a statistical 
one. 

The statistical approach assumes certain averaging. Sometimes, as with disorder, one can think about an ensemble 
of chaotic systems. In such cases ensemble averaging is sufficient. For an individual chaotic system, such as a Rydberg 
atom in a magnetic field, averaging over a wide energy interval is the only choice. 

Quantities investigated in the statistical approach to quantum spectra include various correlators of density of 
states (DoS) v{E) — Tr 5{E — H), where H is the Hamiltonian of the system. Here, it is natural to measure energy 
diff'erences in units of the mean level spacing A. Perhaps the property most frequently studied is the dimensionless 
two-point DoS correlator, 

R2{s) = A^ (u (e + ^a) ^(e-'-/\))- 1, (1) 

where s is the dimensionless energy difference. As mentioned above, for disordered metals the statistical average, 
denoted by (■ • ■) can be performed over different realizations of the random Hamiltonian while, in general, the average 
can be taken over a wide energy band. 

Associated with each particular chaotic system there are typically two relevant energy scales. The first, Ec, is 
associated with the classical time scale Tc = h/Ec on which a density distribution in phase space becomes ergodic, 
i.e. spreads uniformly over the constant energy shell. On time scales larger than Tc, time averages over a classical 
trajectory can be substituted by microcanonical averages over the energy shell in phase space. In a cavity in which 
a quantum particle scatters ballistically from a boundary, the "chaotic billiard" , the energy scale E^ is typically set 
by the frequency of the shortest periodic orbit, or the inverse flight time across the system. In a weakly disordered 
metallic grain, on the other hand, the classical energy scale is set by the inverse transport time, or Thouless energy 
Ec = TiD/L^ , where D denotes the classical diffusion constant, and L represents the system size. The second energy 
scale is set by the mean energy level spacing A which defines the Heisenberg time tr = h/ A. 

The two energy scales can be combined into the dimensionless ratio, 

g = E,/A (2) 

which represents the "dimensionless conductance" ||] of a chaotic system. The ergodic time Tc = h/Ec sets the scale 
beyond which the details of the classical dynamics of the system become irrelevant (as long as the system is chaotic). 
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Correspondingly, for energy scales s ^ g, spectral statistics become universal, independent of the details of the 
underlying classical dynamics. To a very good approximation they are determined only by the global symmetries of the 
system [|[-|| such as T-invariance, and coincide with the universal Wigner-Dyson level statistics of the corresponding 
random matrix ensembles 1^: Unitary (broken T-invariance), Orthogonal (spinless T- invariant systems), or Symplectic 
(T-invariant systems with spin-orbit interaction) |10 In the semiclassical limit, and in dimensions greater than 



one, the dimensionless conductance is large, g 3> 1, and universal statistics apply over a wide energy interval. 

The study of the universal and non-universal statistical properties of quantum chaotic systems has been conducted 
largely along two separate lines: 

1. The first approach has been based on the nonlinear tr- model proposed by Wegner and is applied to the 
study of an ensemble of similar systems such as weakly disordered metallic grains in which electrons experience 
scattering by a random potential. Here, ensemble averaging is a crucial step exploited at an early stage of 
calculation. The supcrsymmetric version of the nonlinear cr-model proposed by Efetov [p^ provided a micro- 
scopic justification for the random matrix theory (RMT) description of universal long-time properties of such 
systems . The same theory accounts for non- universal features of spectral statistics associated with the 
diffusive dynamics of electrons. However, such an approach suffers from two drawbacks: Firstly, it relies on 
the very existence of an ensemble. Very often we are concerned with non-stochastic chaotic systems, such as 
a chaotic billiard, where the notion of an ensemble is inappropriate. Secondly, this type of averaging tends to 
erase information about individual features of the system. 

2. The second approach has been based on Gutzwiller's trace formula in which the semiclassical DoS is expressed 
as a sum over the classical periodic orbits [Q. This approach focuses on the behavior of individual systems, and 
statistical properties reply implicitly on energy averaging. This approach has proved to be a powerful tool in 
describing non-universal properties associated with short-time behaviour. However, its success in reproducing 
universal long time properties associated with RMT has been limited 16|. In particular, it fails to account for 



the correct behavior at times in excess of the inverse level spacing or the Heisenberg time tjj = h/A. Thus, in 
spite of extensive numerical evidence fl^,^ which supports the random matrix description of spectral statistics 
of individual systems at small energies, the origin of its success has remained obscure. 

In Ref. jl8| it was conjectured that the a-model approach to diffusive systems could be generalised to include the 
wider class of chaotic systems. In a recent and insightful development, Muzykantskii and Khmel'nitskii proposed 
a cr-model to describe short-time ballistic dynamics in disordered conductors. Although their argument relied solely 
on impurity averaging, they conjectured that their cr-model should apply even in the absence of disorder. 

Here we develop a semiclassical field theoretic description of quantum spectral statistics of individual chaotic 
systems based solely on energy averaging. This method offers a semiclassical description of the statistical properties 
of individual quantum chaotic systems which accounts both for universal as well as non-universal features. The basic 
ingredients of the underlying classical dynamics are no longer the individual periodic orbits, but general properties 
of the classical flow in phase space. It will be shown the the cr-model constructed in this way indeed coincides with 
that proposed in Ref. |l9| ]. 

Before introducing the main conclusions of this study, we begin by identifying the questions which will be of most 
concern. To do so it is convenient to draw on the insight offered by the study of the dynamics of a particle moving 
in a background of weak randomly distributed scattering impurities. Amenable to the method of ensemble averaging, 
properties of this chaotic quantum mechanical problem are now reasonably well understood. 

What is the quantum evolution of a wavepacket in a background of impurities? According to the time of evolution, 
the dynamics of the wavepacket is characterised by several quite distinct regimes. On time scales t in excess of 
the mean free scattering time r, the initial ballistic evolution of the wavepacket becomes diffusive. At longer times 
the interference of different semiclassical paths induces a quantum renormalization of the bare diffusion constant 
D = VpT/d. This leads to the phenomena of "weak localisation" and is responsible for the quantum coherence effects 
observed in transport properties of mesoscopic metallic conductors. If the impurity potential is not strong enough to 
localise the wavepacket altogether, the wavepacket continues to spread. After a time Tc = LF' / D, the typical transport 
or diffusion time, the wavepacket is spread approximately uniformly throughout the system. Further evolution of the 
wavepacket is therefore said to be ergodic. Beyond the ergodic time t ^ the evolution of the wavepacket becomes 
universal, independent of the individual features of the system. Finally, the spectral rigidity characteristic of quantum 
chaotic systems leads to an approximately coherent superposition or "echo" |^ of the wavepacket at t — th after 
which the wavepacket relaxes to a uniform distribution. 

A similar question can be asked about the quantum evolution of a wavepacket introduced into, say, an irregular 
cavity (quantum billiard) without impurities. In such systems it is widely believed that there too exists some ergodic 
time Tc after which properties of the system become universal. However, at shorter time scales, how is the unstable 
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nature of the classical dynamics reflected in the quantum evolution? Is there an analogue of quantum renormalization? 
What, in general, plays the role of the diffusion operator in describing the low energy degrees of freedom? 

In a recent study by three of us ||2^, a comparison of results taken from the leading order of diagrammatic 
perturbation theory for disordered metals within the diagonal approximation of periodic orbit theory led to the 
conjecture that, for general chaotic systems, the role of the diffusion operator is, in general, played by the generator £ 
of the classical evolution (or Perron- Frobenius) operator, e~^*. If p(x, 0) is an initial probability density distribution 
defined as a function of phase space variables x = (q, p), where q and p are coordinates of the position and momentum 
vectors respectively, the density at a later time t is given by 

p(x, t) = e-^V(x, 0) = ^ - "*(y)] p(y' o), (s) 

where u*(y) is the solution of classical equations of motion with initial conditions y at i = 0, and F covers the region 
of available phase space. One eigenvalue of the Perron-Frobenius operator is always unity and corresponds to the 
ergodic state. The rest of the eigenvalues are of the form e"'''^'*, where 5i(7^) > 0, and appear as complex conjugate 
pairs. They are associated with decaying modes in which a smooth distribution relaxes into the uniform ergodic state 
and are analogous to the diffusion modes of disordered systems. This analogy, therefore, suggests that the ergodic 
time in general chaotic systems is set by the first non-zero eigenvalue Tc ~ ?i/5R(7i). 

As well as confirming the conjecture made in Ref. ||2l|, the field theoretic approach will reveal that, in the semiclas- 
sical limit, all statistical spectral properties of the quantum system depend only on properties of the Perron-Frobenius 
operator. In particular, by considering only its leading eigenvalue 70 = 0, one exactly reproduces RMT, while taking 
into account higher modes enables one to characterize deviations from universality. The field theoretic approach 
provides for the first time a systematic and controlled way to investigate quantum corrections which lie beyond the 
diagonal approximation typically employed in the periodic orbit theory. 

The range of energy scales over which this approach is valid is illustrated in Fig. 1. Energy averaging will be 
performed over a wide energy band of width W centered at Eq. As well as requiring that Eq ^ W, we will assume 
that W is much larger than the energy scale set by the first non-zero eigenvalue of £, 71. This is to ensure that the 
time scale h/W is fine enough to resolve the behavior of the classical dynamics over a time interval smaller than the 
ergodic time ^ fi/ji- It will be also assumed that the finest energy scale, the mean level spacing A, is much smaller 
than 7i. Therefore, we will focus on a range of energy scales where A <C 71 <C <C £^o- 

As well as describing the universal regime of RMT, corresponding to energy scales comparable to the mean level 
spacing A, the field theory developed here properly describes non-universal behavior which appears at larger energy 
scales. However, to avoid encountering non-universal features associated with the finite band width W, it will be 
always assumed that all correlators involve energy differences much smaller than W. Final results are therefore 
expressed in zeroth order in A/W , i.e. when the number of levels in the band tends to infinity. 

The paper is organized as follows: A derivation of the ballistic nonlinear cr-model is presented in section II. The 
interpretation of the resulting functional integral will be discussed in section III. There it is shown that regularization 
of the functional integral forces one to identify the low lying modes of the field theory as the Perron-Frobenius modes 
associated with irreversible classical dynamics. In section IV we present two applications: the reduction of our model 
to RMT, and a calculation of the two-point DoS correlator beyond the universal regime. At the end of this section the 
range of validity of our approach is discussed. In section V we generalize the model to incorporate non-semiclassical 
corrections such as weak scattering from (5-correlated random impurities. This establishes the relation between the 
ballistic cr-model and the conventional diffusive counterpart. A summary and discussion of the results is presented in 
section VI. 
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FIG. 1. A schematic diagram showing the relevant energy scales, as well as the domains of validity of RMT and the field 
theory constructed in this paper. 



II. THE NONLINEAR a-MODEL 



To present the derivation of the effective field theory describing spectral correlations of quantum chaotic systems, 
we will focus on the problem of a single particle confined by an irregular potential described by the Hamiltonian 



(4) 



The classical coimtcrpart of the quantum Hamiltonian is assumc^d to he chaotic and to have no discrete symmetries. 
We confine attention to closed systems so that classical motion inhabits a finite region of the 2(i-dimensional phase 
space. We will assume that all classical orbits are unstable and, in particular, exclude (KAM) systems where the 
phase space contains islands of regular motion. 

We will concentrate on statistical properties defined on an energy band of width W centered at an energy Eq ■ To 
discuss meaningful averages it is necessary to assume that the average DoS, specified by the Weyl formula 



{,^{Eo)) = ^Jd^d[Eo- iJ(x)] , X = (q, p) 



(5) 



is approximately constant within this interval. Taking as an example a particle in a random impurity potential, the 
accuracy of this approximation is of order W/Eq, and can be made arbitrarily small by going into the semiclassical 
limit Eq — > 00. On the other hand, the bandwidth is assumed to be sufficiently large that the number of levels, 
N = h'{Eo)W ^ 1 can be employed as an expansion parameter — final expressions will be expressed in the zcroth 
order approximation in 1/N. Henceforth we will express energy in units of the mean level spacing, A = l/{v(Eo)) 
and denote such energies by e = E/A. For simplicity, it is convenient to employ Gaussian averaging 



(•••)e 



/ 



de 



(27riV2)i/2 



exp 



(•••)• 



(6) 



A general n-point correlator of physical operators, such as the local or global DoS or current densities, can be 

obtained from a generating function which depends on appropriate external sources. Here we focus on two-point 
correlators. Expressed as a field integral, the generating function for two-point correlators takes the form 



{J) = zJ exp J dq*t(q)i (q-i^,) _ fj ^(q) 



(7) 
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where G ^(e) = e — s+A/2 — H denotes the matrix Green function with energy difference s between retarded (R) and 
advanced (A) blocks, J represents the source, and the constant z is included to enforce the correct normalisation.|^ 
Following Ref. j2|] we express the 8-component superfields which appear in Eq. (0), in the block notation 

where superscript p refers to retarded/ advanced components, subscript g refers to fermionic (F) components x 
bosonic (B) components S, and subscript d refers to time-reversal (complex conjugated) components. The introduction 
of equal numbers of bosonic and fermionic fields is a standard trick which obviates the need to introduce replicas and 
normalizes the generating function to unity. Matrices 

A=(^j ®tg®td. ®^d, (10) 

break the symmetry between the advanced/retarded and graded components respectively, and we have chosen a 
convention which introduces the supermatrix 

x=(; i)\u. (11) 

The inclusion of complex conjugated fields effectively doubles the number of fields and implies the relation 

vI/t(q) = *(q)^C^, (12) 
where the operations of complex conjugation, and transposition of supervectors are defined following Efetov p^ , while 

C^t^®(-'I' M , (13) 



ri 



denotes the "charge conjugation" matrix, and 



^^K?o)/ ^^ni-i), 

represent Pauli matrices which act inside time-reversal blocks]^ 

As an application, Eq. (^ can be used to represent the fluctuations in the two-point correlator of DoS by defining 
the source as 

J = J5(q- q')fcA, (15) 
where J is a constant. Then R2{s), defined in Eq. (Ilj), can be expressed as 
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^2(s) = -Y^|^5R(Z(J)),, , (16) 



If the energy difference s is chosen to be much smaller than the width of the energy band TV, correlators become 
independent of N and of the particular shape of the band (whether it is Gaussian or Lorentzian, etc) . Performing the 
energy averaging (o) of Z{J) we obtain 



^ Energy averaging of the generating functional (^) induces a quartic interaction of the form (^'^L^')^ among the supervector 
fields. The matrix L serves as a metric tensor. In Ref. it was shown that the appropriate group of transformations 
preserving the interactions in the fermionic sector is compact while in the bosonic sector it must be chosen non-compact. This 
fixes the definition of L (see Eq. dill)). With this definition, the constant 



z = exp 



isTr, 



In(AL) 



(8) 



accounts for the correct normalisation. For the definition of STrg see discussion below Eq. (|l^. 



This convention differs slightly from that discussed by Ref. |14 



5 



5int = X (/ 



(17a) 
(17b) 

Therefore, in contrast to an impurity averaging, energy averaging induces a nonlocal interaction of 5*. This represents 
an important departure from the usual consideration of random Hamiltonians. 

The next step involves the decoupling of the interaction induced in the averaging by means of a Hubbard- 
Stratonovich transformation. This is achieved with the introduction of 8 x 8-component supermatrix fields (5(qi,q2) 
which are non-local in space. To define the correct decoupling it is crucial to identify those contributions to S'int which 
vary slowly in comparison with the wavelength. In the semi-classical analysis that follows we will show that the low 
lying degrees of freedom are described by matrices (5(qi,q2) which vary slowly with respect to the center-of-mass 
coordinate (qi + q2)/2. Anticipating this, it is convenient to switch to a momentum space representation of the 
interaction and explicitly separate such contributions 



Sint = ^ I dP I dp 



«'t(p)i*(-p)j {^t\p + p)i^'(-p - p 

*t(p)i4»(_p)^ (^*t(_p _ p)L^(p + p 



(18) 



The characteristic momentum cut-off is defined such that po < W/v, where v is the velocity of the particle. Using the 
charge conjugation symmetry of ^' ( p^ ) it is straightforwar d to show that these terms give an equivalent contribution. 
Therefore, following Ref. [Q we decouple the interaction ( |l7b| ) as 



STrq [^ + y*t^Q^ 



(19) 



where STrq denotes the trace operation for supermatrices, STrilf ~ TiMpF — TrMss, with a subscript q used 
to denote a further extension of the trace to include the coordinate integration. Eq. (n2|) implies that the dyadic 
product A(q, q') = ^'(q) ® \I't(q')L obeys the symmetry property A(q, q') — C'^ LA^ (^,q)LC . This induces the 
corresponding symmetry 



Q = C^LQ^LC, 



(20) 



where the transposition should be understood in the sense of an operator. 

Substituting Eq. ( [l9| ) into Eq. (17a) and integrating over ^f* we obtain the following expression for the averaged 
generating functional 



I 



DQ exp 



isTrqO^ + isTrqln (g-\Q) - 



J 



eo-H - iNQ. 



(21a) 
(21b) 



Thus far no approximations 
obtain an effective action. To 
minimizes the action in Eq. (21< 
of Q around the saddle-point, 
massless modes. The integral 
order in 1/A^ (see Appen dix 
symmetry of the action ( [l7al) 
CT-model. 



have been made. The next step is to identify the low energy degrees of freedom and 
do so, we will employ a saddle-point approximation and find the matrix Qo which 
a ) . The effective field theory is described by the expansion of the action in fluctuations 
These fluctuations are strongly anisotropic and can be classified into massive and 
over the former can be evaluated within the saddle-point approximation to leading 
|). The integral over the remaining massless modes, which arise from the underlying 
must be evaluated exactly. The resulting field theory has the form of a nonlinear 



A. Saddle-point approximation and the a-model 



Varying the action in Eq. (21a) with respect to Q, and neglecting the terms of order s and J, we find minima at 
Qo which satisfy the equation 



(22) 
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where Qq must be treated as an operator. The saddle-point solution which is diagonal in superspace is given by 



Qo 



, eo-g 
' 2N 



1 - 



ep-H 
2N 



1/2 



A. 



(23) 



Note that NQq plays the role of the self-energy in the average Green function Q{Qo)- 

The saddle-point solution in Eq. ( p3| ) is not unique but is in fact one member of a degenerate manifold of solutions. 
Their existence follows from the underlying symmetry of the action of Eq. (17a). The interaction term S'i„t['I'] is 
invariant under the group of transformations — > U'^ such that 



U^LU = L, 



(24) 



where U is an operator in Hilbert space. Terms that break the symmetry of the total action in Eq. (17a) are sA, JkA 
and the commutator [H, U]. The invariance of the relation ^E't = '^'^C'^ under the transformation ^E* U'i' induces 
an additional constraint on U 



(25) 



From Eq. (|T^) it follows that these transformations induce the following constraint on the Hubbard-Stratonovich field 

The saddle-point solution in Eq. is not invariant under this group of transformations. Therefore the low energy 
modes of the action are of the form Q = U^^QqU. However, not all of these transformations should be taken into 
account. The group of transformations (|2j) contains a subgroup of matrices which commute with the Hamiltonian. 
The matrix Q remains diagonal in Hilbert space in the basis of eigenstates of the Hamiltonian. In Appendix ^ we 
show that the massive mode integration gives rise to a suppression of the fluctuations of Q by the large parameter N. 
The only matrix U commuting with H which "survives" the N oo limit is the one proportional to the unit matrix 
in Hilbert space. Admitting matrices Q of such form into Eq. (pla|) we obtain 



where 



DQ expi-SeffiQ]), 



(26) 



SeffiQ] = --STrqln 



g-\Q)-—A-J 



g-\Qo) -U {—A + j]u-' - U[H, U-'] 



(27) 



The last three terms under the logarithm in Eq. (27) are small as compared to the first, and an expansion can be 
made in them. Each order in this expansion brings an additional power of 1/N, and suggests the inclusion of just the 
leading order term: 



SeffiQ] = ^ STrq 



Q[^—A + J-U-'[H,U] 



(28) 



This approximation is justified only if s <C and the commutator [H, U] is n ot a nomalously large. The validity of 
this approximation must be considered individually for each system. In section V A we will discuss an example where 
this is not the case, and one has to keep the second order expansion of the logarithm in Eq. (p7|). 



B. Semiclassical approximation 



In the lim it ep — > oo, the configurations of the Q-matrix that contribute substantially to the functional integral 
in Eq. (17a) can be described within the semiclassical approximation. It is therefore convenient to re-express all 
operators in the Wigner representation. Given an operator O as a set of matrix elements C'(qi,q2) between two 
position states at qi and q2 , its Wigner representation is a function of the phase space variables x = (q, p) defined 
by 
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0(x)= J df^q' exp(zp-q7n) 0(q + q72, q - q72) 



(29) 



We will use the fact that, in the semiclassical limit, the Wigner transform of a product of operators is equal the 
product of the Wigner transformed operators, (C'iC'2)(x) Oi(x)C'2(x), where C'i,2(x) are smooth slowly varying 
functions on the quantum scale [E3| . In this approximation Eq. (EJ) becomes 



U\-k)LU{'x) = L, 



(30) 



and implies that the matrices U(x) belong to the pseudounitary supergroup C/(2,2/4). Expressed in the Wigner 
representation, the constraint in Eq. (25), 



f/*(q,p) =CC/(q,-p)C^ 



(31) 



shows that the matrices C/(x) at different x are not independent. This agrees with the findings of Rcf. [p4[ . 

The massless modes in the Wigner representation are generated by those matrices C/(x) which do not commute 
with A. Such matrices, denoted by r(x), belong to the coset space H = G/K = [/(2,2/4)/ [U{2/2) x U{2/2)] and, 
as implied by Eq. (p^), satisfy the symmetry relation 



T*(q,p) =Cr(q,-p)C' 



(32) 



As mentioned above, those matrices U which commute with the Hamiltonian are strongly suppressed by massive modes 
(see Appendix ^). In the semiclassical limit, we therefore admit only those matrices T'(x) which are independent of 
the energy. The massless modes are then given by 



Q(x) = r-i(x||)go(F)T(x||), 



(33) 



where X|| denotes a phase space coordinate on the energy shell eq — H{x.). 

Substituting T for U in Eq. ( |2^ ) , and applying the semiclassical approximation in which the commutator with the 
Hamiltonian becomes the Liouville operator C, 



[H,f] -inCT{^\\) = -i?i{T(x||),ff} , 



where 



dA dB dB dA 

dqi dpi dqi dpi 



(34) 



(35) 



denotes the Poisson bracket of A and _B, we obtain 



Q(x) ( + J + ihT-^CT 



(36) 



Since the only dependence on the coordinate x±_ — -ff(x) normal to the energy shell enters through Qo{H), given by 
Eq. (p3|), the integral over this variable can be performed and yields a factor ttNA. Introducing the notation 



Q(x||) ^^JdH T-i(x||)Qo(i?)r(x||) = r-i(x||)Ar(x||), 
we obtain the final expression 



(37) 



Q ( + J + ihT-^CT 



(38) 



Here and henceforth when the arguments of Q and T are omitted they should be understood as functions of X|| 



Note that, since we adopt the convention in which the DoS is equal to unity, the phase space coordinates are normalized as 

/ dX||//l'' = 1. 
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The matrix Q(x|[) introduced in Eq. (|37|) satisfies the nonUnear and symmetry constraints 



Q(x||)' 



1, Q(q,p) = C^LQ^(q,-p)LC. 



(39) 



Naively, since Eq. ( |38| ) is expressed through the matrices T rather than through Q = T^^AT, it appears that it 
does not correspond to a cr-modeL However, a general property of <7-models is the invariance of the action under 
gauge transformations T RT, where R commutes with A. Using the fact that the Liouvillian £ = X|| • Vxy (where 
X|| is the phase space velocity) is a first order differential operator it is straightforward to show that, under a gauge 
transformation, the change of the action (^) is given by 



6Seff[Q] = 



STr 



KR-^CR 



■nfi 



(ix|| 
17 



STr 



KC Ini? 



= 0. 



(40) 



To arrive at the last equality we used the fact that the flow in phase space is incompressible: Vxy • xy = 0. 

The kinetic part of the action in Eq. (^) is equivalent to that introduced by Muzykantskii and Khmel'nitskii |li 
where it is written in the form of the Wess-Zumino-Witten action 



wzw 



[S] 



7r?i 



du STr Q 



dQ 

du 



(41) 



where Q is a smooth function of the auxiliary variable u and X|| with the boundary conditions: Q(x||, 1) = Q(x||) 
and Q(x||,0) = A. The equivalence of Eq. ( p] ) with the kinetic part of the action in Eq. ( |38| ) can be established 
straightforwardly by substituting Q = T^^AT and manipulating the various terms using the identities T'^CT = 
-{Cf-^)f and f-^duf = -{duf'^)f which follow from C {f-^T^ = du {f-^Tj = 0. An integration by parts with 
respect to X|| shows that the resulting integrand is a total derivative with respect to u. The integration over u, with 
the boundary conditions T(x||,l) — r(x||) and T(x||,0) = 1, leads to the kinetic part of the action in Eq. (|3^). 



C. Range of validity of the a-model 

To clarify the domain of applicability of the nonlinear cr-model in Eq. (|3q) let us review the main steps involved in its 



derivation. The construction of the effective generating functional in Eq. (21a) involved purely formal manipulations 
which involved no approximation. To proceed beyond this expression we invoked a saddle-point approximation in 
which the fluctuations of the massive modes were neglected. The parameter which controlled this approximation was 
the inverse bandwidth 1/A^. 

The second approximation involved the replacement of quantum mechanical commutators by the semiclassical 
Poisson bracket. Such an approximation is justified at high energies where the shortest length scale is set by the 
wavelength of the particle. Finally, in treating fluctuations of the massless modes around the saddle-point, we treat 
only the leading order term in the expansion. Formally, if the commutator [_ff , T] is not anomalously large, this 
approximation is also justifled by large N . Since characteristic configurations of T are assumed semiclassical this 
assumption can be violated only if H contains some non-semiclassical contributions. 



The validity of this semiclassical approximation is discussed in more detail later in section V A when we return to 
consider scattering from quantum impurities and the relation of the ballistic cr-model to the conventional diffusive 
nonlinear cr-model. 

The symmetry breaking terms in the action place additional constraints on the range of validity. The expansion 
around the saddle-point relies on characteristic frequencies (or energy scales arising from the Poisson bracket in the 
action (^sj)) being much smaller than the bandwidth N . 

The derivation of the cr-model of Eq. ( ^8|) relies solely on the presence of energy averaging which allows us to neglect 
the contribution of massive modes in the functional integral. Indeed, energy averaging is crucial in the ballistic limit 
even in the presence of disorder. This was emphasized in the study of Altland and Gefen of spectral statistics of 
ballistic metallic grains. There it was pointed out that ensemble averages of spectral correlators differ from averages 
performed over both ensemble and energy. In the semiclassical language of periodic orbit theory this difference emerges 
from trajectories which are not scattered by impurities [p6| and give rise to "clean features" in the quantum spectrum. 
The neglect of interference terms among different trajectories (namely the diagonal approximation) is allowed only 
upon energy averaging over a wide band. Otherwise, the interference among these trajectories is substantial. In 
section V A we show that without energy averaging only the diffusive cr-model of Efetov can strictly be justified. 



In this case, the large parameter which suppresses the fluctuations of the massive modes is (rA) 
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As numerical or experimental studies of quantum chaos always involve finite statistics it is worth discussing how 
the results of such studies should be compared to the expressions above. Since the exact DoS consists of i5-functions 
at the positions of energy levels, all spectral correlators based on finite statistics are inevitably singular. As the 
statistics are increased by extending the energy window W over which averaging is performed, these correlators 
remain singular. The meaningful way to define spectral correlations is by introducing a smoothed DoS, p-yi^e) — 
(7/77) l/((e — e^)^ + 7^). For any finite 7, the limit N ^ 00 exists and generates a "smooth" function. As an 
example, consider the two-point correlator of DoS for the unitary random matrix ensemble: To obtain the universal 
expression i?2(s) = 6{s) — sm^{TTs)/ {ns)^ it is possible to take the limit 7^0 only after the limit — > 00. If, 
however, one keeps N finite while decreasing 7, the expression will approach the universal result only until 7 « 1/A^ 
after which it will start to deviate with increasing magnitude. The expressions for spectral correlators that we are 
describing here should be understood in exactly the same way. They are asymptotic expressions corresponding to the 
limit of A — !■ cx) taken before the limit 7 — > 0. When one deals with experimental/numerical data, which necessarily 
involve finite statistics, one should always keep the level width finite such that 7 > 1/A^. 



III. REGULARIZATION OF THE FUNCTIONAL INTEGRAL AND IRREVERSIBILITY OF THE 

CLASSICAL DYNAMICS 



In this section we discuss the question of ultraviolet divergences of the cr-model and show that the regularization 
procedure forces us to understand C in Eq. (pq) as t he g enerator of irreversible classical evolution. 



The functional integral in Eq. (|26| ) with the action (38) suffers from ultraviolet divergences and must be regularized. 
The ultraviolet divergence is not an artifact of the approximations employed in the derivation of Eq. ( ^8[ ) and is present 
even in the original expression for the ratio of quantum spectral determinants Det(e— -ff). Although supersymmetry 
of Eq. improves the ultraviolet properties of the functional integral, in higher dimensions it is not sufficient to 
make the expression converge. Therefore an ultraviolet regulator needs to be introduced. This regularization induces 
the corresponding regularization on the functional integral in Eq. (p6|). The kinetic part of the effective action ( |3^ ) 
can be presented as STr[QT~^X|| VxyT], where X|| is the classical phase space velocity. This action is only sensitive to 
the variations of the Q-matrix along the classical trajectories. Therefore nothing prevents the Q-field from fluctuating 
in the directions transverse to X|| . It is these short scale fluctuations that ultimately lead to the divergence of the 
functional integral. The ultraviolet divergences are independent of the classical dynamics and of the shape of the 
constant energy surface and are unphysical. The diverging contribution to the functional integral therefore needs to 
be extracted by an appropriate regularization procedure. 

The problem of the ultraviolet regularization of functional integrals is well studied in field theory. One of the ways 
to regularize the functional integral is to introduce a term 

6SR = mJ dx||STr(Vx|| Qf, (42) 

into the effective action (|3^). This term suppresses strong fluctuations of Q in the directions transverse to X||, and 
favors the physical functions Q(x||) which are smooth. Depending on the dimensionality of the phase space this 
may not be sufficient to make the integral convergent and additional regularization procedures should be invoked. To 
explore this issue in more detail we will consider the functional integral which arises from the lowest order perturbative 
expansion of the action. 

In this case we can represent T = 1 + ST and expand the action (^8|) with the regulator (^) to second order in 
ST. We refer the reader to section IV B| where this is discussed in greater detail. Here we only outline the conceptual 



steps which relate to the regularization. 

In the lowest order of perturbation theory, the resulting Gaussian functional integral generates simply the deter- 



minant of the elliptic operator is — Cr = is — C — mW'^^^ (see section IV B). An operator is called elliptic is the 
component of highest rank in derivatives is positive definite. The problem of regularization of the determinants of 
such operators is discussed in the literature (see, for example, Ref. [0). One method involves the construction of a 
zeta function of the operator defined as 

1 1 r / \ 1 

C(^. - Cn\z) = E (— ^ - ^ / .-^-dm-exp [- - £,) t] . (43) 

Here denote the eigenvalues of Cr, and we assume that s is chosen such that the operator is — Cr has no zero 
modes. Then the integral in the right hand side converges at the upper limit. At the lower limit i ^ it can diverge 
depending on the value of z. However, this divergence is ultraviolet in nature and can be removed by taking the 
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integral at sufficiently large positive z. The expression can then be analytically continued to the rest of the complex 
plane. 

A regularized spectral determinant is expressed through the derivative of the zeta function 

lnDet(is-£fl) = -C'(i5-^j?|z) . (44) 

The regularized spectral determinant Det(«s — L^) is a function of s. It has zeroes at the positions of eigenvalues of 
Lb. and nowhere else on the complex plane. By taking the limit one recovers the result which is independent of the 
regulator (|2|) 

— ^ = Det(is - t)^ lim Det(is - Lb)- (45) 

Z\IS] m^O 

This limit is very different for integrable and chaotic systems. In particular, antihermiticity of L in the limit m — > 
suggests that the zeroes of the regularized determinant Det(is — C) lie on the imaginary axis of is. However, for 
chaotic systems this is not the case |p^ , p9| . 

To understand the subtleties which arise when this limit is taken for nonintegrable systems, let us consider the 
purely classical evolution. Suppose we form an initially non-uniform probability density distribution p(x||) in the 
phase space. The classical dynamics involves stretching along the unstable manifold and contraction along the stable 
one. Thus, any non-uniform initial distribution will evolve into a highly singular function along the stable manifold. 
The regularization term (^) in the classical evolution can be ignored for short times but eventually, when contractions 
along the stable manifold make the phase space gradients sufficiently large, it becomes relevant.^ Therefore the limits 
time-to-infinity and m ^ do not commute. To find the spectrum one has to take the time-to-infinity limit first and 
then set the regulator to zero. In this limit the eigenvalues 7^ of L have finite real parts corresponding to relaxation 
rates into the equilibrium distribution. These physical eigenvalues which reflect intrinsic irreversible properties of the 
•purely classical dynamics are known as Ruelle resonances or the Perron- Frobenius spectrum [ p8| , p9[ . 

Thus, the necessity to regularize the functional integral of the tr-model forces us to understand C as the classical 
evolution operator which corresponds to the irreversible classical dynamics. 

There are several ways of calculating the Perron- Frobenius spectrum other than diagonalizing and taking the 
"zero noise limit", m — ^ Oj3l|. These employ, for instance, symbolic dynamics |]3^ , course graining of the flow 



dynamics in phase space and analytic continuation [ p4[ . An exact formal expression for the dynamical zeta 
function 1/Z{z), which should be understood as a regularized product Jl^^C-^ ~ 7m)' given in terms of the classical 
periodic orbits of the system. For two-dimensional systems it is of the form BE 



where Tp is the period of the p-th primitive orbit and A is the eigenvalue of the Monodromy matrix with absolute 
value larger than one. (The Monodromy matrix is the linearized map on the Poincare surface of section in the 
vicinity of the orbit.) In its present form, 1/Z{z) cannot be used to determine the eigenvalues 7^;. For this purpose 
a re-summed formula is required. It can be obtained by expanding the infinite product over the periodic orbits and 
ordering the various terms in a way that leads to maximal cancellation among them. This method, known as the 
cycle expansion p6[ | , exploits the property that the dynamics of chaotic systems in phase space is coded by a skeleton 
of a small number of periodic orbits called fundamental orbits. In this sense, long periodic orbits may be viewed as 
linear combinations of the fundamental orbits. 

To summarise the main conclusion of this section, it was shown that the regularization procedure which properly 
defines the functional integration forces one to understand the low lying degrees of freedom of the action ( p8|) as the 
Perron- Frobenius relaxational modes of the classical counterpart. We emphasize that these modes represent purely 
classical characteristics of the system independent of the regularization procedure. Indeed the corresponding spectral 
determinant det(2: — C) has an exact representation in terms of the classical periodic orbits of the system (Mq). 



*An analogous situation arises in the theory of turbulence [^. In the inertial range, viscosity can be neglected and turbulence 
can be considered as dissipationless. However at sufficiently small scales, velocity gradients become large and viscosity becomes 
relevant. In this picture the energy which is pumped into the system at large spatial scales is transfered without dissipation in 
the inertial range to smaller spatial scales and is eventually absorbed at microscopic scales determined by viscosity. The latter 
can be viewed as an ultraviolet regulator, which is eventually set to zero but has a finite effect on the velocity correlators since 
it is necessary to produce a stationary solution. 
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IV. APPLICATIONS 



To interpret the findings of the previous sections we will apply the generalized nonlinear cr-model to the regime of 
long-time or low energy scales. This will establish a firm connection of level statistics with RMT. Corrections to RMT 
will be studied within the framework of a perturbation theory involving the modes of the Perron-Frobenius operator. 
These results indicate a close correspondence between spectral correlations of the classical and quantum operators 
which we discuss. 



A. Random Matrix Theory 

It is widely believed that the statistical quantum properties of systems with few degrees of freedom can be described, 
at least over some range of energy scales, by RMT 17pj7|. To interpret this, various approaches have been developed 
largely along two parallel lines discussed in the introduction. The first approach concerned the study of ensembles 
of random systems such as disordered metallic grains [^,^ 38 1. Randomness in this case is introduced on the level 



of the Hamiltonian itself usually as a consequence of some impurity configuration. The second approach involves the 
study of non-stochastic systems which are chaotic in their classical limit such as the Sinai or the stadium billiards . 
In this case "randomness" is generated by the underlying deterministic classical dynamics itself. Nevertheless, it has 
been conjectured |^] that spectral fluctuations of strongly chaotic quantum systems are described by level statistics 
of random matrix ensembles. 

Despite being supported by extensive numerical studies, the origin of the success of RMT as well as its domain of 
validity are still not completely resolved. Below we will show that, in the semiclassical limit, this conjecture is indeed 
valid for chaotic systems without any discrete symmetries, and which are characterized by an exponential decay of 
classical correlation functions in time. 

If we define by {7,1} the set of eigenvalues of the Perron-Frobenius operator £, then the lowest eigenvalue in ergodic 
systems is 70 — 0. This eigenvalue, associated with the invariant density on the energy shell, is non-degenerate 
and manifests the conservation of probability density. Any initial density distribution eventually relaxes to the state 
associated with 79. If, in addition, this relaxation is exponential in time, then the Perron-Frobenius spectrum has a 
gap associated with the slowest decay rate. Thus, for the first non-zero eigenvalue 71, we have 7^ = 5R(7i) > 0. This 
gap sets the ergodic time scale, = I/71 over which the classical dynamics relaxes to equilibrium. In the case of 
disordered metallic grains, it coincides with the Thouless time, while in ballistic systems or billiards it is of order of 
the time of flight across the system. 

In the limit s <C 7^, or equivalently at times which are much longer than t^, the dominant contribution to Eq. ( |2^ ) 



with the effective action of Eq. (38) arises from the ergodic classical distribution, the zero- mode CTq — 0. Taking only 
this contribution, the functional integral ( p6| ) becomes definite. Previous studies have demonstrated the equivalence 
of the zero-mode action with Wigner-Dyson level statistics of RMT |l4|JT5|] . For example, taking the source as given 
by Eq. (|l5|), we obtain 



Z{J)}eo = J dQo exp [-i JSTr ( [s+ + 2Jk] AQo)] , (47) 



where Qq = T^^ATq. This expression coincides with that obtained in Ref. [|T^ and reproduces Wigner-Dyson level 
correlations. We therefore conclude that the quantum statistics of chaotic systems having no discrete symmetries and 
with exponential classical relaxation are described by RMT at energies smaller than 7J . 

The RMT description is expected to hold even for certain chaotic systems where the Perron-Frobenius spectrum is 
gapless [ 3|. Examples include the stadium or Sinai billiards where classical correlation functions decay algebraically 
in time |3^. In this case, the resolvent {z — jC)~^ is expected to have cuts which reach the axis. Nevertheless, we 
expect the RMT description to hold whenever the spectral weight of the resolvent inside the strip < Sftz < 1 (which, 
however, excludes the pole at the origin) is much smaller than unity. 



B. The two-point correlation function: Beyond universahty 

In this section we will make use of the cr-model to examine how corrections to RMT appear at larger energy scales. 
Again, focusing on the two-point DoS correlator, the generating function leads to the expression 



i?2(s) = j^^JvQ (^y"dx||STr[AfcQ(x,[)]j exp[-^eif(s)] 



(48) 
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where 



Seffis) 



(ix||STr 



i—AQ + QT-^CT 



(49) 



Although straightforward, the perturbative expansion is somewhat technical, and here we present only the results of 
the detailed calculation described in Appendix ^ 

In the limit of high frequencies s 3> 1 , the two-point correlator takes the asymptotic form 



R2{s)^Rp{s)+Rnp{s), 



(50) 



where both the nonperturbative term Rnp{s) as well as the perturbative one Rp{s) are expressed through the classical 
spectral determinant 2?(s) as 



RNpis) 



cos(27rs) 9, , , , ^ ^ , , 

—^V\s), Rp{s) = --—\n[V{s 



(51) 



The determinant 2?(s), regularized according to the procedure outlined in section II], is expressed in terms of deter- 
minants of the Perron- Frobenius operator 



P(s) = 5ft- 



Det'(£)2 



Det 



{is — £)(— is — C) 



(52) 



where the prime indicates that the zero eigenvalue should be excluded from the determinant. 'D(s) can be expressed 
in terms of the eigenvalues 7^ of £, the Ruelle resonances, as 



hi + s') 



2 ' 



(53) 



where ^(7^) = 7^ for 7;^ 7^ and ^(70 = 0) = 1. Note that, if the product in Eq. ( [s^ ) is formally divergent, I?(s) 
should be understood as the regularized determinant (^ . 

These results agree with those co nject ured in Ref. | pi| | and compare with the perturbative expressions previously 
found for weakly disordered metals |p8| , p^pc| ] when the eigenvalues of the Liouville operator are identified by the 
eigenvalues of the diffusion operator. 



V. BEYOND THE SEMICLASSICAL APPROXIMATION 

The derivation of the nonlinear cr-model in Eq. ( ^8|) relied on the use of the semiclassical approximation. However, 
often we are concerned with quantum chaotic systems which can not be treated straightforwardly within the framework 
of semiclassics. A familiar example involves the quantum mechanical scattering of particles from a weak random 
impurity potential. In such cases, a formal justification of the ballistic nonlinear cr-model in Eq. ( |3^ ) does not seem 
possible. However, if the quantum Hamiltonian can be resolved into a part that can be treated within semiclassics 
and a part which can not, when the latter is small, a perturbation treatment may still be possible. 

Consider a general Hamiltonian H with matrix elements 

H = H,i + Hqu, (54) 

where Hd represents the contribution which can be treated within a semiclassical approximation, and i/qu determines 
the part which can not. 

If the matrix elements of Hq^ are small as compared to the band width N (a more precise criterion for a specific 
operator Hqu is formulated below) their effect can be treated within the tr-model approach. In this case the saddle- 
point is governed by Hc\, and we can use Eq. (|2^) with H replaced by Hd- The contribution of Hq^ to the effective 
action can be found by expanding Eq. (|27|), 



SeffiQ] = -^STrqln 



g-'{Qo) -Hq^-U[—A + J]U-'- U[H, [/-ii 



(55) 
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where the supermatrix Green function involves only Hd- Expanding to second order in Hq^ we obtain 

1 ' ■ ' 



SeffiQ] = 2^ STrq 



2N 



(56) 



Finally, representing the Q-matrices in the Wigner representation, the second order correction to the action takes the 
form 



where i?qu(q,q') = (q|iJqu|q'}. 



3_ f lTrf„ TT ^^p-'Pl(qi-q2)/2?i-ip2(q3-q4)/2?l 

4iV2 J 11^9^1 /jd^ 

^ — 1 ^ — 1 



j=l i=l 

xi?q„(q2, q3)i/qu(q4, qi)STr [Q (pi, (qi + q2)/2) Q (pa, (qa + q4)/2)] , 



(57) 



A. Random impurities and the restoration of the diffusive nonlinear a-model 



To illustrate these ideas, let us consider the physical example involving a particle moving in a background of weakly 
scattering impurities. If the Q matrices vary on a scale that is long as compared to the scattering length £ = vt, 
the particle dynamics becomes diffusive and we should recover the supersymmetric nonlinear cr-model obtained by 
Efetov [p^ . In the opposite limit, the impurities generate a new term in the action which takes the form of a collision 
integral. 

The problem of dilute scattering impurities in an otherwise ballistic system has been discussed previously. A 
description within the framework of diagrammatic perturbation theory was investigated by Altland and Gefen [ p5| . 
More recently, in an important development, Muzykantskii and Khmelnitskii introduced an effective field theory 
to extend the diffusive a-model into the ballistic regime. 

For simplicity, let us consider a (dimensionless) (5-correlated white noise Gaussian random potential 



i?qu(q,q') = ^(q)'5'(q-q'), 
with a second moment defined by the mean free time r, 

h 



(<5F(q))^ = 0, {6Vi<DSVic^))y 



(58) 



(59) 



Here (• ■ ■)v denotes the ensemble average over the random potential, u = 1/Ail represents the average local DoS, 
and Q is the volume of the system. 

In this case, the expansion of the action around the saddle-point of the Hamiltonian Hd is justified in the limit 
h/r ^ NA. The same condition allows the truncation of the perturbation series at second order. Once again, 
performing the energy integration ( |3^ ) and using the fact that dx\\/h'^ — dqdp^^/Anp'^fl, where py is momentum on 
the constant energy shell (|p||| —pp), we obtain the effective action 



dqdpil 



Str 



dqdpiidp'i 











-A + J 


-ihT 









{inp%y 



STr 



Q(q,Pii)Q(q,p;) 



(60) 



Although this action is precisely of the form of that introduced in Ref . |l9[| , its derivation and the domain of validity 
seems far removed from that proposed in this earlier work. The cr-model description of the ballistic regime holds only 
if the frequencies of interest (or, equivalently the characteristic energies of the gradient terms) are small as compared 
with the width of the band iVA. In the absence of energy averaging the range of validity of this description is restricted 
to the diffusive regime, where it coincides with the diffusive cr-model [Q. At higher energies the massive modes have 
to be taken into account, and the cr-model description breaks down. The distinction drawn by energy averaging has 
been emphasized by Altland and Gefen |p5| . Physically the difference arises from those orbits whose period is shorter 
than T but longer than the inverse band width (A^A)"^ Technically the energy averaging suppresses the massive 
mode fluctuations and facilitates the cr-model description. 

To establish the relation between the ballistic cr-model and the conventional diffusive counterpart we follow Ref. . 
Let us suppose that the classical component of the Hamiltonian corresponds to free propagation. Anticipating a rapid 
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relaxation of the momentum dependent degrees of freedom of Q on the energy shell, and a slow variation of the spatial 
modes, we introduce a parametrisation which involves the moment expansion 



T(x) = rA'(q)To(q), TK{(i) - exp[in(q) • K(q)], 



(61) 



where n = p/|p| and, without any loss of generality, we choose AK + KA = 0. To enforce the nonlinear constraint in 
the most convenient way, we have adopted a parametrisation which departs from that discussed in Ref . ]l9| . 
Expanding the action to second order in K and performing integrals over n we obtain 



'S'cff = 



dqSTr 



iA 



A + J]Q- . {Vn)T,'A Ak 



where Q(q) = Tq ^(q)ATo(q). Performing the Gaussian integration over K we obtain the effective action 



riD{VQf + i2A ( s+A + 2 J ) Q 
which coincides with that of the conventional diffusive cr-model [[14 



(62) 



(63) 



VI. DISCUSSION 



In conclusion, we have shown that the quantum statistical properties of chaotic systems are described by a functional 
supersymmetric nonlinear ct- model with an effective action given by Eq. (38). This result was obtained by employing 
energy averaging as opposed to ensemble averaging previously used for disordered metallic grains. The low lying 
degrees of freedom of the action ( |3^ ) were identified as the Perron-Frobenius eigenmodes of the underlying classical 
dynamics. Thus, statistical characteristics of the quantum mechanical system in the semiclassical limit are determined 
by the symmetries of the system and properties of the Perron-Frobenius operator. In particular, a universal behaviour 
described by RMT is expected whenever the system has no symmetries and a gap exists in the Perron-Frobenius 
spectrum. 

Our approach, however, assumes no systematic degeneracies of the actions of the classical orbits of the system other 
than those associated with the known discrete symmetries of the Hamiltonian. To emphasize this point, consider 
the return probability to a given point. In the semiclassical limit it is given by a double sum over classical returning 
trajectories Y^^ - MA*, where Ai denotes the probability amplitude associated with the z-th path. If the corresponding 
actions are much larger than ?i, the return probability reduces to a sum over probabilities, rj^^ l^iP- The factor 
rj is an integer which accounts for exact degeneracies in the orbits of the action, which arise, for example, from the 
existence of time reversal or reflection symmetries. When such a degeneracy is characterized by the existence of 
a discrete symmetry, it can, in principle, be incorporated into the cr-model. However, there are systems in which 
degeneracies of the actions cannot be characterized by simple discrete symmetries. Examples include the arithmetic 
billiards on surfaces of constant negative curvature [Q. The actions of the periodic orbits of these systems become 
exponentially degenerate as their length increases. This is a result of hidden symmetries which originate from number 
theoretic properties of these billiards. Indeed, despite showing a gap in the corresponding Perron-Frobenius spectrum, 
arithmetic billiards to not exhibit random matrix behaviour. 

The cr-model derived here has a wide domain of validity which goes well beyond the RMT results. Using perturbation 
theory the two-point DoS correlation function (|^) was calculated. The result (shown in Eqs. (^0|) and (pi])) is expressed 
in terms of the spectral determinants of the Perron-Frobenius spectrum. Similar results were obtained recently by 
Bogomolny and Keating However, their results differ from ours by terms which are related to high repetitions 
of the same periodic orbit. The two results therefore clearly coincide in the limit where all orbits are highly unstable. 
At this stage the source of discrepancy is not understood. Whether it is related to corrections to the leading order 
of perturbation theory (Eqs. (^0|) and (^l|)) or to the nature of uncontrolled approximations used by Bogomolny and 
Keatin g ||42| remains unclear. It is appropriate, however, to mention that the cr-model functional integral with the 
action (|38|) can be solved exactly for the one-dimensional harmonic oscillator [ 0| and gives the correct result. This 
result can also be obtained by use of the perturbation theory described in section IVB. Such a calculation shows that 
it is necessary to take into account the contributions of both stationary points of the action to obtain the exact result. 
The diagonal approximation of periodic orbit theory, on the other hand, also reproduces the exact result. Since all the 
periodic orbits of this system are repetitions of one primitive orbit, this result suggests that the diagonal approximation 
commonly employed in periodic orbit theory, does not coincide precisely with the diagrammatic perturbation theory, 
as is commonly assumed. 
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Many of the results previously obtained for disordered systems concerning statistical properties of wavefunctions 
and spectra depend only on the spectral properties of the diffusion operator in a given system Q,^. These can be 
generalized straightforwardly to the case of chaotic systems by substituting the spectrum of the diffusion operator by 
the Perron-Frobenius spectrum. 

The field theoretic approach described in this paper offers a systematic way of studying a variety of issues. These 
include; (i) the transition between the orthogonal and the unitary ensembles in ballistic systems; (ii) the effects 
of discrete symmetries on spectral statistics in systems exhibiting hard chaos; and (iii) weak localization effects in 
ballistic systems p6| . 
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APPENDIX A: SADDLE-POINT APPROXIMATION: IDENTIFYING THE MASSIVE MODES 



In this section we examine fluctuations around the solution (^3|) of the saddle-point Eq. (^2|) to identify the massive 
modes in the effective theory of Eq. ( |21b|) . Although we focus our remarks on the orthogonal case studied in this 
paper, the general conclusions of this appendix hold for all ensembles. 

To identify the massive modes it is convenient to work in the eigenbasis {<p„} of the quantum Hamiltonian, where 
Eq. (p3|) takes the form 



2N 



1 - 



ep - 
2N 



1/2 



A 



(Al) 



Massive modes appear as fluctuations SQ that commute with Qq in superspace. Expanding the action in Eq. ( pib| ) 
around the saddle-point to second order in SQ and neglecting s (and using the fact that Qq is diagonal in the Hilbert 
space indices) we obtain 



(A2) 



The mass of these modes is not apparently large but is of order one. However, their contribution to the two-point 
DoS correlator is given by 



R 



1 



2, massive 



(47riV) 
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1 



(A3) 



where (• ■ •)g denotes the average over supermatrices Q with respect to the action in Eq. (21a). This contribution 
vanishes in the N ^ oo limit. 

If we consider the contribution of the massive modes to a correlator of local observables such as the local DoS i^(q) 
we find 



1 



massive 



^2' ^ 



¥'M(qi)¥':(qi)^;(q2)^.(q2) (STr(Afc5g^,) STr(Afc<5Q,^))( 



(AttN) 

■^^'pMi)^*^^) X ^(^*(qi)(/3^(q2) 



(A4) 



The contribution from both the diagonal = v) and off-diagonal terms is small: The former is of order N^^ , while 
the latter involves N"^ terms each of which is of order N^"^ . However, since the the phases of wave functions at different 
point are almost uncorrelated so the off-diagonal terms arise with random phases. This implies a contribution of the 
off-diagonal terms which is also of order N^^. Another way to see this is by invoking the completeness argument: 
Each sum in the last line of Eq. (A4) tends to 5{c\i — q2) as the width of the band is increased (this follows from 
completeness of the basis of eigenstates of the Hamiltonian). At any finite band width the function (y3^(qi)iy9* (q2) 

has the characteristic width [VL/NY^'^, where £7 is the volume of the system, and can be approximated by the Heaviside 
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function {N/n)Q{{Q/Ny^'^ — |qi — q2|)- Therefore the r.h.s. of Eq. (A4) vanishes if the coordinates qi and q2 are 
separated by a distance larger than {fl/NY^'^. These considerations enable us to neglect the massive modes. 



The integration measure in Eq. ( |21a| ) is invariant under the group of transformations Q U ^QU, where U is an 



operator satisfying Eq. ( p4| ) with indices both in the Hilbert space U^i, and in superspace. The action Eq. (21b) is 
also invariant under such transformations, provided that U commutes with H. We will denote such transformations 
by Uo This symmetry leads to the existence of a degenerate manifold of saddle-point solutions (at s — 0). All matrices 
of the form 

Q - U^'QoUo, (A5) 
where [IJo,H] = satisfy Eq. (22). In the basis of the eigenstates of the Hamiltonian such matrices are of the form 



Uo.fj.1, — Sfj_„Uo,fj_ with C/o,M S UOSP{2, 2/4) |15|. We assume the absence of degeneracies due to non-integrabihty. All 
such matrices generate zero-modes. 

It is shown below that the integration over the massive modes strongly favors the ground state configurations of 
Q which correspond to identical Q^^'s. This happens because the ground state in which Q^'s are different breaks 
supersymmetry of the action for the massive modes. This leads to a rapid decay (as a function of inhomogeneity of 
Q^) of the superdeterminant which arises from the integration over the massive modes. Therefore the integration 
over the massive modes gives a non- vanishing contribution to the effective action which depends on Q^'s. This 
contribution can be interpreted as an effective interaction between Q^'s which favors configurations with identical 
Q^'s. Hence, it can be thought of as "ferromagnetic" interaction of "spins" which reside on the nonlinear manifold 
C/OS'P(2,2/4)/[[/OS'P(2/2) x UOSP{2/2)]. This interaction is long range (all "spins" within the band interact with 
approximately equal strength) and therefore, in the thermodynamic limit N oo, leads to a ferromagnetic ground 
state. The fluctuations of "spins" from the the ground state configurations are small as 1/iV and can be neglected. 

To see how the supersymmetry breaking for the massive modes arises let us consider one term in the sum ( [A2| ) 
corresponding to particular ij, and u. The matrix has the same symmetries as the Q-matrix in Efetov's nonlinear 
(T-model and can be parametrized as 



-isin(6'p)) -cos(6'^) 

Here we deviate from Efetov's original parametrization by introducing the angles 6*+ and 6~ which can be expressed 



through the angles appearing in Ref. |14|, 6i and 62 as 9^ = 9i + 92 and 9~ = 9i ~ 92- The particular form of the 
matrices u and v is not important for what follows and will be left unspecified. 

If the angles 9^ and 9^, coincide then the massive modes line up with 9^. In other words we can make a global 
rotation to bring 9i, to zero, and in this coordinate frame the massive fluctuations correspond to SQ^^ and SQ^j^^. If 
the angles 0^ and 9^ differ by a small amount, we can go to the "center of mass" coordinate where 

= (A7) 

In this frame the massive modes will still correspond to SQ^^ and SQ^^. The contribution of 5Q^^ to the effective 
action is 

STr {QE{9t,W5Q^^u,QE{9.)u,5Qt^u^ + 5Q^^6Qt^) (A8) 

Instead of integration variables 5Q^^ and SQ^^ it is more convenient to use 5Qf^^ ~ u^SQ^^Uy and SQ^^ = 
UySQ^^u^. Since the superjacobian of such transformation is equal to unity. 



SDet = 1. m 

the invariant measure is preserved. 

With the parametrization involving ordinary variables a^, 6^, and Grassmann variables Ui, ct*. 



/ ai 02 icji i(T2 \ / «i —0.2 —0-3 -cta 



AA 




6Qt^=\ I, (AlO) 



CTg CT4 ih\ ib2 I ' I —icrl ~i(J2 ib\ 162 

V (74 CTs ib2 ib\ I \ —iu\ —iai 162 ibi 
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which obey the symmetry relations 



integration over massive modes 5Q^^ and SQ"^^ with the action Eq. (A8) can be performed and yields 

[2 + cos6l(cosh6l+ + cosh6i-)]4 



'-fiiy 



(All) 



(A12) 



(2 + cosh^ e+ + cosh^ e-){2 + 2 cosh 6*+ cosh 6*-) (2 + 2 cos2 0)^ ' 
For small 6*^, writing cos0 = 1 — a, cosh6'+ = 1 + and cosh6'^ = 1 + (3^ , /^j/ can be expanded to second order, 



l--(2a + ^++^-)2 



exp 



1 



-(2a + /3+ +/3^)2 



(A13) 



If all 9^ are small, then we obtain a model equivalent to spins with infinite range interactions. In the thermodynamic 
limit of such a model the mean field approximation becomes exact. The fluctuations of a, (3^ and /3 ~ be come small 
as 1/A^ and can be neglected. This forces us to consider the matrices Qo which are of the form of Eq. (Al). Then the 
relevant (massless) fluctuations of the Q-matrix are those that anticommute with A in superspace. 



APPENDIX B: PERTURBATION THEORY 



In this appendix we will employ the cr-model to study to the two-point correlator of DoS fluctuations. In particu- 
lar, we will examine the perturbative corrections to RMT which appear at larger energy scales. To obtain the high 
frequency asymptotics of i?2(s) the functional integral in Eq. (^) can be evaluated using the stationary point approxi- 
mation. However, to obtain the contribution which is non-perturbative in 1/s it is necessary to introduce an additional 
term u^{AQ)'^ into the action (^9|) which serves as a regulator controlling the stationary point approximation |p^ , ^ . 
Ultimately, the regularization parameter u can be set to zero. 

We therefore express the two-point correlator as 



i?2(.s) = lim 

li— >o 54 



VQl / dx||STr[AA;g(x||)] exp[-S'ej(s 



where 



dx||STr 



i—AQ - u^AQ f + QT-^CT 



(Bl) 



(B2) 



The derivation of the results already presented in Eqs. (^ ) and (|5^) closely parallels that of Ref. |l^. For a more 
detailed account of the method see Ref. . At high frequency the integrand in Eq. ( p^ ) becomes highly oscillatory, 
and we can use the stationary phase method to evaluate the integral. We will show that there are two stationary 
points: Q = A and Q = —Ak. The term u^{AQ)'^ in the action is introduced in order to stabilize the second one. We 
can expand the integrand in small fluctuations of the Q-matrix around A and — Afc to obtain the leading high-frequency 
asymptotics of i?2(s)- 

We first consider the expansion around Q = A. This corresponds to the ordinary perturbation expansion previously 
employed in the study of disordered conductors |l4,pq]. We begin with the parametrisation 



T = l + iP, 



P = 



q B 

B 



where, from Eq. (|3^), it follows that P satisfies the condition 

P(q,p)* = -CP(q,-p)C^. 



(B3) 



(B4) 



Next we substitute Eq. (B3) into Eq. (Bl) and expand the integrals in the pre-exponential factor and the free 
energy (B2) to second order in P. Due to the presence of the infinitesimal imaginary part in s^, the stationary point 
Q = A is stable and we can safely set it = in the free energy (B2). To second order in B and B we have 



STr(AfcQ) 
STr(AQ) 



8 - 2STr(fcSB + kBB), STr[AkQf « ~%STr{kBkB + BB) 
-4STr(BB). 



(B5a) 
(B5b) 
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Using these relations we obtain the following expression for the perturbative part of i?2(s) 



Rp{s) = / V[B,B] ( / dx[|[l - -STiikBB + kBB)] ) cxp[-S'e^(s) 



where 



Seffis) = in / (dx||) STr 



-sBB - iBCB 



0{B^ 



In order to perform the integration over B and B it is convenient to represent these matrices as 

3 3 



(B6) 



(B7) 



(B8) 



i=0 



i=0 



where the Pauli matrices are defined in Eq. (|14|). As foUows from Eq. ( pO[ ) the matrices B and i? in Eq. (B3) obey 
the relation B = kB\ which implies 



B, = kB], i = 0, ...,3. 



In this notation, Eq. (B7) becomes 

Seff,p{s) = -ITT j (dx||) STr 
Each matrix Bi can be parametrized as 



(B9) 



(BIO) 



(Bll) 



The parametrization for Bi can be obtained from Eq. (B9). To evaluate the integral (B(;) over the variables (Bll) 
one can use Wick's theorem. It is necessary to take into account Eq. (B4) which reduces the number of independent 
integration variables by a factor of two. As a result we obtain the second part in Eq. (pi]). 

The stationary point Q = A of the functional integral (Bl) is not the only one. To find the other stationary 
points consider Eq. (Bl). It is possible to parameterize fluctuations around a general stationary point Qa as Q = 
)o(l + «^o)(l — iPo)^^ where Pq anticommutes with Qo and no longer obeys equations (B3) and (B4). Expanding 



the effective action in Eq. (B2) in powers of Pq we would obtain the stationarity condition dSejf{s)/dPo ~ 0. 

This route however is inconvenient since the parametrization of Pq will depend explicitly on Qq. Instead it is 
convenient to perform a global coordinate transformation on H = U{2, 2/4)/[J7(2/2) x U{2/2)], Q — > Uq^QUq, where 
C/q £ H, which maps Qo to A. 

Since all points on a symmetric space are equivalent by definition, this coordinate transformation preserves the 
invariant measure and leaves the functional integral in Eq. (Bl) invariant. The integrand, however, will change because 
it contains matrices A and —kA that break the symmetry in the coset space. Such a coordinate transformation is 
equivalent to changing only the source matrices A Qa = UqAUq^ and — fcA ~QkA ~ — C/q/cAC/q"^ in Eqs. (pl|), 
(B2) and keeping the old parametrization of Eq. (B3). The stationary points will correspond to those Uq for which 



the linear in P terms in the expansion of the effective action vanish. 

Note that because the transformation matrix Uq is independent of momenta p, from Eq. (^Tj) it follows that Uq — 
CUqC'^. Therefore Uq belongs to the coset space of the usual orthogonal ensemble U0SP{2, 2/A)/[UOSP{2/2) x 
UOSP{2/2)]. For this case there is only one other stationary point ||l^,|^ corresponding to Qa = UqMJ^'^ — —kA 
and —QkA = —UokAU^^ = A. For this point Eq. (Bl) can be rewritten as 



lim — 5 

M^o 64 









1 exp 






CT 





(B3) 



We now expand Eq. (B12a) in powers of P using Eq 

around —kA in Eq. (Bl). E xpan ding the free energy (B12b) to second order in P we use Eq. (B5a) 
the parametrization of Eq. (Bll) 



(B12a) 
(B12b) 



This expansion is equivalent to expanding the Q-matrix 

Note that with 
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STr {kBB + kBB) = -4^(|a,p - \b,\^), 

i=0 
3 

STr {kBkB + BB) = -4^(|a,p + 



(B13a) 
(B13b) 



Th erefore the Grassmann variables in the parametrization (Bll) do not couple to s and u . As follows from 
Eq. (B13a), the ordinary variables and hi couple to s with opposite sign. Due to the prese nce of the infinitesimal 
imaginary part in s the integral over (the zero mode variable) would diverge at u = 0. Eq. ( Bf 3b ) shows that the 
term STr(fcAQ)^ makes the integration over a? convergent. We therefore have to keep u finite during the evaluation 
of the functional integral and take the limit m only in the final expressions. The quadratic approximation to the 
free energy (Bf2b) becomes 



eff[ 



-2TTis + 2TT / (dxji)STr 



■ 3 

E 

.i=0 



C)bi + cr*Cai + r]*Cr], 



(B14) 



The zero mode Grassm ann variables ry^ and do not appear in the quadratic expansion of the effective action (Bf4). 
For the integral ( Bf 2a ) not to vanish they have to come from the pre-exponential factor. While evaluating the integral 
we have to take into account the symmetry (B4) which reduces the number of independent integration variables by 
factor of two. Therefore there are eight independent Grassmann variables in the zero mode. Thus, in order to obtain a 
non-zero result we should expand the pre-exponential factor to eighth order in P. Then in the eighth order expansion 
of the prefactor we should keep only the zero mode terms. This renders the integration over the zero mode variables 
non- vanishing, whereas the integration over the ordinary zero mode v ariab les yields a factor (s^ -t-u"*)"^. The integral 
over the non-zero modes yields the superdeterminant of the operator (Bf4). After we perform the integration we take 
the 7i — > limit to obtain Eq. (^ 
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